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Abstract Generalized Ray Theory (GRT) provides a simple description of MHD 
mode transmission and conversion between magnetoacoustic fast and slow waves 
and is directly applicable to solar active regions. Here it is tested in a simple two- 
dimensional, isothermal, gravitationally-stratified model with inclined magnetic field 
using previously published exact solutions and found to perform very well. 
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1. Introduction 



Magnetohydrodynamic (MHD) mode conversion has long been proffered as an ex- 



planation of observations (Braun, Duvall, and LaBonte, 1988 Braun, 1995) of /- 
and p-mode "absorption" by sunspots. Several methods have been used to explore 
this question, including numerical solution of a differential eigenvalue problem {e.g., 



Cally, Bogdan, and Zweib"ell [19941 [SHgdan and Call^ [19971 |Crouch and Call^[2003 



2005 



2003 



1997 



2008 



with results directly compared to the solar data in [Gaily, Crouch, and Braun 
andlGrouch et al. 20051, and direct numerical simulation (Gaily and Bogdan, 



|Ga ly, 2000 1 jKhomenko and GoUados, 2006 [ [Gameron, Gizon, and Duvall, 



However, a very different technique, based on a generalization of ray theory ( Gaily, 



2006 Schunker and Gaily, 2006) presents several novel insights which greatly enhance 



interpretation and understanding. To date, this has been applied only in two dimen- 
sions (2D), in the sense that gravity, the magnetic field, and the direction of wave 
propagation all lie in the same plane, and we shall maintain that restriction here. 
This decouples the Alfven wave from the problem, leaving only the fast and slow 
magnetoacoustic waves (although see Gaily and Goossens 2008 for a quantitative 
estimation of 3D coupling). 

In particular. Generalized Ray Theory (GRT) has verified that mode transmis- 
sion/conversion occurs at or near the equipartition level (2;oq) where the sound and 
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Alfven speeds coincide, c = a. A magnetoacoustic wave incident on this level is 
partially transmitted (transmission coefficient T) and partially converted (coeffi- 
cient C). GRT introduces the concept of the attack angle (a), the angle between 
the wavevector and the magnetic field at the mode conversion level. If a. is small, 
transmission dominates, but if it is large, then conversion is favoured. 

In any plausible sunspot model, the sound speed (c) varies slowly with height (z) 
through the surface layers, but the Alfven speed (a = B / ^JJip) increases rapidly due 
to the decreasing density (p). Here, /i is the magnetic permeability (47r x 10^^ in SI 
units). For z -ti Zcq we have c ^ a, whilst a ^ c in the opposite regime z 3> ^oq- 
In c ^ a the fast wave is predominantly acoustic in nature, and the slow wave is 
largely magnetic. On the other hand, where a ^ c the opposite pertains: the fast 
wave is magnetic and the slow wave is acoustic. In these asymptotic regimes, the 
fast and slow magnetoacoustic waves are completely decoupled and distinguishable. 
However, in the conversion region (a w c) they interact. Importantly, the "fast" and 
"slow" waves do not maintain their identities through this coupling region. 

Let us be clear about the meaning of "transmission" and "conversion" here. We 
say that a wave is totally transmitted (T = 1) if it maintains its acoustic or magnetic 
identity across the conversion layer. For example, a vertically propagating sound 
wave in vertical magnetic field does not interact with the magnetic field, since it 
is longitudinal, and stays a sound wave as it passes through Zcq. we call this total 
transmission. Conversely, an incident acoustic wave in c ^ a which passes through 
to become purely a magnetic wave is said to have been totally converted. Of course, 
in general, both transmission and conversion are partial. 

GRT is only an approximate description of wave behaviour. It shares many of 
the fiaws of standard MHD ray theory (Weinberg, 1962 I, though not the gross error 
of presuming perfect fast-fast or slow-slow connectivity across Zeq (i.e., T = 0). 
Specifically, it should be asymptotically correct in the high frequency limit, where 
wavelengths are vanishingly small compared to background inhomogeneity length 
scales. But how well does it perform in more realistic, moderate-frequency scenarios? 
In this paper, we test the accuracy of the GRT estimation of transmission coefficient 
T in a simple, uniform-field, isothermal, gravitationally-stratified model for which 
exact wave solutions exist. Although not exhibiting all of the features that we may 
wish in a sunspot model, it does possess the most important characteristic - a rapidly 
increasing a/c ratio with height - which allows for a meaningful and informative test. 



2. Model and Equations 

We adopt the simplest model which exhibits the required features of increasing a/c 
with height, with a = c at some finite level: an isothermal, gravitationally-stratified 
atmosphere with uniform magnetic field inclined at angle 9 to the vertical. Without 
loss of generality, the equipartion level (zcq) is set at zero. A sound wave with horizon- 
tal and temporal dependence exp[i(fca; x — ujt)] is launched upward from great depth, 
is incident on 2; = where it is partially transmitted and partially converted, and 
the resulting transmitted sound wave is examined in 2 > 0. The aim is to calculate 
the transmission coefficient (T) at varying frequencies (lu), horizontal wavenumbers 
{kx), and field inclinations (0). In order that the sound wave is vertically propagating 
rather than evanescent, and that it is indeed an acoustic and not a gravity wave, 
{iu,kx) is chosen to lie in Region I of the acoustic-gravity wave dispersion diagram 
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Table 1. The a-coefficients which appear in equation ([sjl for the exact transmission 
coefficient T. 
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(Cally, 2001 Figure 1), i.e., above the upper branch of uP' —Lo^—(?kl.{\—N'^ /lo^) = 0, 
where oJc = c/2H is the acoustic cut-off frequency, N'^ = g/H — jf? is the square 
of the Brunt- Vaisala frequency, g is the gravitational acceleration, and H is the 
density scale height. This requires u) > bOc for kx = 0, and a slightly more stringent 
condition for nonzero kx- However, in strong, inclined magnetic field, the "ramp 



effect" reduces the effective cut-off frequency to oJc cos & ( 


Schwartz, Cally, and Bel, 


1984 


Schunker and Cally, 2006 


Mcintosh and Jefferies, 2006 


1. This effective cut-off 



2.1. Exact Solution 



Following |Cally| ( |200T| [2008t , we define the following dimensionless variables: fre- 
quency V = ujH/c, vertical position C, = ujH/a where H is the (constant) density 
scale height, and the Brunt- Vaisala frequency n = NH/c = (7 — l)^^^/7, where 7 
is the usual ratio of specific heats. In these units, the acoustic cut-off frequency is ^. 



Since a = ce^/'^" it follows that C = ve ^^^^ and that C 
K = kx H is, the dimensionless horizontal wavenumber, kq 



•0+ 



as 2 — 

2 oOi^2 1 



00. Finally, 
- 1/4 , and 

Kz = \/ v'^ -\- {n^ — v'^^K^jv'^ — 1/4 is the vertical dimensionless wavenumber in the 
equivalent non-magnetic case. The definition of kq directly invokes the ramp effect, 
as it is real only for u) uj^ cos^ Q. 

With these variables, the linearized MHD wave equations may be expressed as 
a single fourth order ordinary differential equation of hypergeometric type for u, 
the component of velocity perpendicular to the field. The resulting exact magneto- 
acoustic-gravity wave solutions may be expressed in terms of Meijer-G functions 



(Zhugzhda and Dzhalilov, 1984a I, or more simply the 2^3 hypergeometric function 



(Cally, 2008 1. The 2^3 function is extensively discussed in Luke (1975). 

In the degenerate case of horizontal magnetic field, the equation reduces to second 
order, and the well-known spectral structure discussed at length in Chapter 7 of 



Goedbloed and Poedts ( 2004 ) results, including cusp (often called slow) and Alfven 
continua. These continua are due to the coefficient of the highest (i.e., second) deriva- 

and Goedbloed 



five vanishing at Alfven and cusp "critical levels" (see Cally 



1984 



and Poedts, 2004 Section 7.3.2). There are no such critical levels when the magnetic 



field has a vertical component, and they will not be discussed further. 
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As set out by Cally (20081, the solution for u is 



2 1^ ~\~ i/^2 : 



1 — 2k, ^ — k — ikq — Ik tan 0, ^ — k + ikq — ik tan 0; — ("^ sec^ 9^ + 

C2 C^" 2-F3(^ + K - iHz, ^ + K + IKz; 

1 + 2k, i + k — Iko — ZK tan 6, i + k + Ikq — Ik tan 0; — sec^ + 



1 — 2iKo+2iH: tan i 



2-^3(1 — Ikq — Ikz + iKtan^, 1 — Ikq + ikz + Ik tan ^ 



1 — 2iKo, f — Iko — k + Ik tan 0, | — Ikq + k + Ik tan 0; — sec^ 0j + 

C4 ^l+2i'^0+2iK tan 9 ^^g^l -(- J/^q _ i/j^^ ^ IK tan 0, 1 + Ikq + IKz + iKtanfl; 

1 + 2iKo, I + Ikq — k + iKtanff, | + Ikq + k + iKtan^; — ("^ sec^ 0) , 

(1) 

where the C,; are the arbitrary constant coefficients. The individual solutions ui, 
. . . , Ui multiplied by the C-coefficients are said to be of Types 1 to 4 respectively. 
Since the 2^3 functions all approach 1 as — > (z ^ 00), the asymptotic behaviour 
of these four solutions is evident, and they may be identified as the exponentially 
growing fast mode, the evanescent fast mode, the outgoing slow (acoustic) mode, 
and the incoming slow mode respectively. Applying top regularity and radiation 
boundary conditions, the first and fourth solutions are therefore dropped, so the 
physical solutions are necessarily a linear combination of U2 and U3. Similarly, we 
may apply a bottom boundary condition [z — > —00, i.e., C, — > 00) that there be 
no incoming slow modes (allowing only an incoming fast mode and outgoing slow 
and fast waves). |Cally| ( |2008| shows that the physical solution for this case of an 
incident acoustic (fast) wave from below is 023^2 — 022^3, where the a-coefficients 
are expressed in terms of F-functions (Table [2?T| ). 

The total wave energy flux and direction can be established for each wave type. 



with the vertical flux of each wave mode remaining constant with height ( Zhugzhda 



and Dzhalilov, 1984a I. The incident fast wave is purely acoustic at large height, so 



neglecting the magnetic term the vertical flux is 



Fz = Re[pi v*z] = Re 



-i K^3 



+ 2^^ 



(2) 



where pi is the Eulerian pressure perturbation, Vz is the z-component of the plasma 
velocity, and and ^z are deflned as the horizontal and vertical components of 
displacement respectively. 

The exact transmission coefficient (T) is defined as the proportion of incident 
wave energy flux transmitted from fast to slow acoustic waves. It takes a similar 



form to that found in the vertical field case by Cally ( 2001 1, and with the appropriate 
a-coefficients is 



T - 



\a22\ (f) 



|a2 3a4 2 - 022043!^ / 



(3) 
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Here the incident fast flux is 



f = b2i K Kz 



4(7 — VjK^v^ + 4ri^ ^1 + 7(n^7 — 1 — 71^^)^ 



+ 



- (7 - l)Ki/{-An-K - + + - 4i/^) 



sin 20 



cos^ e + A{-i -1)k,zV^{i/ - K^)sin0 I , (4) 



where 



62 



f '2. '2. 22 4 22 4 

(7 — l)(riK — 2kj^ + +(riK — J^) cos 26 



+ 71^(7 — 2)^K^i/^ sin^ 20 



Similarly the vertical transmitted flux is 



63 "I 7^1/* + ^1 +7(27;^^ ~ "'^•^) (^n^ + 1^^) cos 26 — 2kkqv'^ sin 20^ 
48ko - 3 + 64(^2k'' + 2j/'* + /«^(8i/^ - 1)^ - 8 ^8k^(1 + 4i/^) cos 261 
3i^^ sec^ + Qv'^ sec"* 6 + 4k«o (^1 + 8k^ + 4ko + (1 + 4ko) cos 26^ sin 26*^ | , 



X KqV 



(5) 



( '2, { '2, '2, 2 4 22 4 \^ 2242 

64 7 I n K — 2kj^ + + {n k — v ) cos 20 1 + (7 — 2) k sin 20 



where 
1 

63 ^ 



2.2. Ray Solution 

Standard ray theory treats waves like a particle moving in phase space. It is applicable 
in weakly inhomogeneous media in the high frequency limit. The 2D magneto- 
acoustic-gravity wave dispersion function on which it may be based is 

T> = T>2D = uj^ — {a^ + c^)fc^a;^ + a^c^fc^fc^ + (? N'^k^ — (lo^ — a^k^ cos^ 6)uj^ , (6) 

where k = {k^ + k^)^^^ = y is the magnitude of the wavevector, kx and kz are 
its horizontal and vertical components, and fcy its component along the magnetic 
field direction. Locally, the dispersion relation D = restricts the solutions within 



frequency- wavevector phase space to fast and slow hypersurfaces (Schunker and 
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Cally, 2006 I. The ray paths are derived from the Hamiltonian equations 



dx 

d7 



dV 



dk 

d7 



dV 
dx. ' 



dt 
d7 



dV 



duj 
d7 



dV 



(7) 



where r is a time-like parameter following the disturbance. In the present case, the 
dispersion relation is independent of horizontal position, (x), and time (t), meaning 



that frequency {uj) and horizontal wavenumber (kx) are constant (Weinberg, 1962). 
The inclusion of acoustic cut-off and Brunt- Vaisala frequencies means that the ray 



paths no longer obey Fermat's principle (Barnes and Cally, 2001). 



Using the foundations set out by Weinberg ( 1962 1 and Tracy, Kaufman, and 



Brizard (2003), Cally (2006) and Schunker and Cally (2006) developed CRT to allow 



for mode conversion within the ray description. Standard ray theory breaks down 
at close avoided crossings of the fast and slow phase loci in z-kz phase space, but 
CRT redresses this by effectively doing a wave-mechanical matching across these 
regions. This process then returns the (approximately) correct connectivity. Readers 



are referred to Schunker and Cally ( ,2006 ) for a full description. We simply precis the 
mathematical process here, without further elaborating on its foundations. 

The mode transmission/conversion capabilities of CRT are based on constructing 
the dispersion function 



P = det D = DaDb 



from a Hermitian dispersion matrix 



D 



Da V 



(8) 



(9) 



Specific expressions for Da and Dt are given in Schunker and Cally (2006). Now, 77 is 



"small" away from the conversion/transmission region, and so the dispersion relation 
= D « Da Dh yields Da = or Df, = 0, the decoupled acoustic and magnetic 
waves. However, close to the "star point" where Da = -Dh = 0, the coupling term 
77 dominates, and provides the local connectivity description. The fraction of energy 
transmitted from acoustic-to-acoustic or magnetic-to-magnetic is given by 



T = Tf = exp 



-27r|7y|7|B| 



where 



B={Da,Db} 



dDa dDh dDb dDa 



dkz dz 



dkz dz 



(10) 



(11) 



is the Poisson bracket of the uncoupled dispersion functions. The subscripted star 
indicates that the expression should be evaluated at the star point where Da and 
D}j simultaneously vanish. Hence I??!* = —T)*. Clearly, T = 1 if 77* = 0, representing 
total transmission. The subscript "f" in Tf is used here to distinguish the "full" form 
of the transmission coefficient from the "simplified" form, which we now introduce. 

Figure [l] displays the avoided crossings in z-kz phases space for typical vertical 
and inclined field cases, and clearly illustrates how magnetic field inclination can lead 
to much higher (and lower) transmission coefficients by narrowing (or widening) the 
gap between the fast and slow loci. 
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Figure 1. Phase space diagrams (z-k^) for the case re = 0.5, v = 1.5, 7 = 5/3. Left: 9 = 0°; 
right: d = 25°. The fast and slow loci (full curves) are labelled in the left panel. The Da = 
and Dj, = loci are represented as dashed curves. The star points are where these cross. In 
the left panel Tf = 0.63 at both star points (the vertical field case is always symmetric), whilst 
in the right panel Tf = 0.93 for the close avoided crossing and Tf = 0.19 for the wide one. 



At high frequency, uic and N can be neglected, moving the star point exactly to 
a = c, and yielding a simplified and intuitively useful expression for the transmission: 



exp 



|fc.|(fc2 + A;2) 



(12) 



where h = [d{a?' / (?) / dz\a=c is the equipartition layer scale height. Note the depen- 
dence on the perpendicular component of the wavevector fc^ = fcsina: when the 
wavevector is parallel to the field it vanishes, and transmission becomes total in the 
simplified description. This is not exactly the case for Tf. In the current isothermal 



context. Equation (12 1 takes the form 



exp 



IT u \K cos ( 



(13) 



The full transmission formula for Tf retains the acoustic cut-off and Brunt- Vaisala 
frequencies, and so should be more accurate, although also more complex. Both Tf 
and Ts are tested against the exact transmission T in Section [3] 



3. Results 

Figure [2] displays T, Tf, and Ts as functions of dimensionless horizontal wavenumber 
(k) for magnetic field inclined at 40° to the vertical and various frequencies u ranging 
from only slightly above the acoustic cut-off at = ^ and Brunt- Vaisala freqtiency 
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Figure 3. Exact (solid line) and full GRT (dashed line) maximum transmission points plotted 
against i/ hi = 40° (left) and the magnetic field inclination angle for i/ = 1 (right). 



n = 0.4899, where 7 = | is used throughout, up to = 10 = 20i^c- As k increases, 
so does the inchnation of the wavevector; at a particular vahie it is paraUel to the 
magnetic field, resulting in k± = k± H vanishing and Ts peaking at 1. The attack 
angle (a) is also displayed on the top axis. (The upper limit attained by k in these 
graphs is due to the model leaving Region I of the acoustic-gravity propagation 
diagram beyond that point.) 

As expected, both approximate formulae perform very well at high frequency. The 
simple ray formula of course incorrectly returns a peak transmission of 1 at a = 0, 
but remarkably the full ray formula manages to almost perfectly fit the exact peak, 
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Figure 4. Exact (solid line), full GRT (dashed line), simple CRT (dotted line) transmission 
formula plotted against k with u = 1 for 9 = (top left), 6 = 20° (bottom left), 6 = 40° 
(top right) and 9 = 60° (bottom right). Accuracy of both GRT formulas decreases with high 
magnetic field inclinations. 



even at very low frequency. Figure |3] shows just how well the peak transmission is 
approximated by Tf for a range of frequencies and magnetic field inclinations. It 
is not surprising that even Tf fails at low frequency and large attack angle, as the 
foundational transmission formula ( 10 1 is derived from a local analysis around the 
star point that becomes progressively more inaccurate as the avoided crossing widens 
(attack angle increases). 

Figure |4] similarly fixes frequency at = 1 and plots T, Tf, and Ts against k (or 
a) for various field inclinations. Both Tf and Ts cope with near-vertical field better 
than highly inclined field, although the peak is well-fitted by the full ray formula 
throughout. 

As an aside, we mention that the exact transmission curves presented here are 
not consistent with Figure 2 of Zhugzhda and Dzhalilov ( 1984b I, where transmission 
coefficient peaks are typically much lower. However, the beautiful correspondence 
between the completely independent exact and full GRT solutions here provides 
confirmation of these present results. 



4. Conclusion 

Generalized Ray Theory (GRT) redresses one of the major failings of classical ray 
theory in the MHD context; the failure to recognize wave transmission/conversion 
at avoided crossings of the fast and slow loci in phase space. It does this by adding a 
local wave-mechanical analysis in the neighbourhood of these "star points" in a way 
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that is quite convenient for practical calculation. The applications to helioseismology 
by Schunker and Cally (2006 I and Cally (20071 have provided important new insights 
into the physics of near-surface mode mixing, which complement lessons learned from 
direct numerical solution of the wave equations (e.g., Cally| 2000 Khomenko and 



CoUados 20061. 



However, it is important to realize the strengths and limitations of this new tool. 
As with all ray theory, it works best at small wavelength. However, we have seen here 
that it is nonetheless extremely useful even at frequencies comparable to Uc and A''. 
The nature of the isothermal model precluded our reducing the frequency further, 
below Uc, as the required incident sound wave then does not propagate at all. 

It is important to appreciate though that the true value of CRT lies not in its 
ability to yield accurate transmission coefficients, which it certainly can do near the 
transmission peaks, but rather in its qualitative and physically compelling description 
of the processes in play. 
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